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INTRODUCTION 

Pressure  fluctuations  have  been  observed  to  occur  in 
most  fluidized  beds  and  these  fluctuations  have  been  used 
to  define  an  index  for  the  quality  of  fluldizatlon  (Shuster 
and  Kisliak,  1952;  Sutherland,  1964;  Winter,  1968).  Small 
and  rapid  fluctuations  are  considered  to  be  associated  with 
a  good  quality  of  fluldizatlon.  It  appears  that  analyzing 
the  relationship  between  pressure  fluctuations  and  bubble 
behavior  in  a  gas-solid  fluidized  bed  gives  rise  to 
effective  information  on  control  strategies  for  the  system 
(Jones  and  Pyle,  1971;  Whitehead  et  al . ,  1977).  The  nature 
of  pressure  fluctuations  in  a  fluidized  bed  is  a  complex 
function  of  particle  properties,  bed  geometry,  pressure  in 
the  bed,  and  properties  and  flow  conditions  of  the 
fluidizing  fluid. 

Pressure  fluctuations  have  been  studied  by  numerous 
investigators  (Tamarln,  1964;  Swinehart,  1966;  Hiby,  1967; 
Kang  et  al. ,  1967;  Llrag  and  Littman,  1971;  Wong  and  Balrd, 
1971;  Fan  et  al. ,  1981).  They  have  drawn  essentially  two 
different  conclusions  about  the  cause  of  pressure 
fluctuations  in  a  fluidized  bed.  Tamarin  (1964)  and  Hiby 
(1967)  have  concluded  that  the  pressure  fluctuations  are 
related  to  the  passage  of  bubbles  through  the  upper 
boundary  of  the  bed  and  to  the  changes  in  the  height  of  the 
bed.  Similarly  Llrag  and  Littman  (1971)  have  revealed  that 
the   fluctuating  signals   of  pressure  involve  the  periodic 
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component  resulting  from  the  bubbles  escaping  from  the 
surface  of  the  fluidized  bed.  Kang  et  al .  (1967),  however, 
have  concluded  that  the  action  of  bubbles  causes  changes  in 
the  mode  and  condition  of  gas  flow  and  porosity  in  the 
dense  phase,  inducing  the  pressure  fluctuations.  Hlby 
(1967)  also  has  proposed  a  deterministic  model,  which 
recognizes  the  presence  of  periodic  components;  the 
frequency  of  the  whole  bed  is  estimated  from  the  assumed 
frequencies  of  individual  particles.  Moritomi  et  al. 
(1980)  have  established  a  deterministic  model  taking  into 
account  inertia,  damping  and  frictional  forces. 

Statistical  methods  were  employed  fairly  extensively 
by  previous  investigators  (see,  e.g.  Swlnehart,  1966;  Kang 
et  al . ,  1967;  Lirag  and  Littman,  1971;  Pan  et  al,  1981) 
used  statistical  approaches.  Swlnehart  (1966)  calculated 
the  cross-correlation  function  between  two  pressure 
fluctuation  signals  taken  from  two  vertically  separated 
pressure  taps  and  determined  the  "correlation-average 
propagation  velocity"  of  bubbles  in  a  fluidized  bed.  Kang 
et  al.  (1967)  and  Lirag  and  Littman  (1971)  calculated  the 
probability  density,  auto-correlation  function,  and  power 
spectral  density  functions  of  the  pressure  fluctuations 
off-line.  The  resultant  auto-correlation  and  power 
spectral  density  functions  were  then  used  to  determine  the 
frequency  of  the  fluctuations.  Lirag  and  Littman  (1971) 
also  used   the  measured   frequency  to  estimate  the  bubble 
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size.  These  investigators  have  amply  demonstrated  the 
power  of  the  statistical  analysis  or  correlation  approach 
in  analyzing  the  random  pressure  fluctuation  signals  in  a 
fluidized  bed;  nevertheless  relatively  little  has  been  done 
to  stochastically  model  the  phenomena  of  pressure 
fluctuations. 

Hiby  (1967)  has  predicted  theoretically  that  frequency 
should  decrease  with  increased  velocity  but  has  failed  to 
prove  it  experimentally.  Jones  and  Pyle  (1971)  found  that 
frequency  increased  with  the  superficial  velocity,  u,  to 
the  undamped  natural  frequency  of  the  system;  but,  no 
effect  of  u  was  found  for  particles  with  a  size  of  500  fi. 
Lirag  and  Littman  (1971)  observed  that  the  frequency 
Increased  slightly  with  the  superficial  velocity  except  in 
the  shallow  beds  of  500  n  glass  particles.  The 
experimental  or  theoretical  results  of  the  various 
investigators  appear  to  contradict  each  other  either.  This 
contradiction  need  be  resolved,  which  is  the  objective  of 
the  present  study. 

In  the  work  reported  in  this  paper,  pressure 
fluctuations  due  to  the  bubble  movement  near  the 
distributor  plate  were  measured  by  means  of  a  pressure 
probe  over  a  range  of  time.  The  resultant  time-series  have 
been  analyzed  off-line  by  determining  their  auto- 
correlation function  and  power  spectrum  with  emphasis  on 
the  effect   of   the  superficial  velocity  of  the  fluldizlng 
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gas.  A  stochastic  model  of  bubble  motion  In  a  fluldlzed 
bed  has  been  developed.  This  model  visualizes  the  bubble 
motion  in  a  fluldlzed  bed  to  consist  of  the  random 
movement,  generating  irregular  signals,  and  the  linear 
movement,  generating  wave-like  signals.  A  theoretical 
auto-correlation  function  and  the  corresponding  power 
spectral  density  function  have  been  derived  based  on  the 
model;  these  functions  are  compared  with  the  experimental 
results. 


THEORETICAL 

In  a  gas-solids  f luidized  bed,  bubbles  appear  at  any 
locale  in  the  bed.  The  pressure  fluctuations  caused  by 
bubble  movements  can  be  measured  by  means  of  a  pressure 
probe  (Lirag  and  Littman,  1971;  Fan  et  al. ,  1981);  the 
detailed  bubbling  phenomena  can  be  characterized  by  the 
auto-correlation  function  and  the  power  spectrum  of  the 
signals  from  the  probe.  Moreover,  stochastic 
characteristics  of  the  bubble  behavior  can  be  reasonably 
determined  from  analyzing  these  auto-correlation  functions. 
The  auto-correlation  function  is  a  manifestation  of  the 
change  in  the  void  fraction  in  the  bed  and  contains  two 
basic  components;  one  essentially  is  a  regular  periodic  or 
wave- like  component  generated  by  the  passage  of  bubbles 
through  the  probe  field,  and  the  other  the  random 
component,  due  mainly  to  the  stochastic  nature  of  bubble 
motion  into  and  out  of  the  field.  The  random  component 
distorts  the  wave-like  component.  The  power  spectrum 
yields  Information  related  to  the  local  velocity 
distribution  of  bubbles  and  the  major  bubbling  frequency. 

Each  recorded  signal  is  assumed  to  comprise  two 
components.  The  first  is  the  periodic  or  sine  wave 
component,  Y(t),  measured  from  its  mean  value.  This  can  be 
expressed  as  (Yutani  et  al. ,  1983) 


Y(t)  ■  A  sin(2jrfot  +  6)  (1) 

where 

A  *  amplitude 

f   ■  frequency  in  cycles  per  unit  time 

e  ■  initial  phase  angle  with  respect  to  the  time 
origin  in  radians 
Since  a  bubble  enters  the  field  of  a  probe  in  a  random 
manner,  6  is  a  random  variable;  it  is  assumed  that  e  is 
uniformly  distributed  with  a  probability  density  function 
of  p(e)  over  (0,2ir). 

The  second  component  is  represented  by  a  stochastic 
process  describing  the  random  passage  of  bubbles  through 
the  field,  for  which  we  define 


lo, 


1,  if  a  bubble  is  in  the  field 
Z(t)  -  1  (2) 

otherwise 


Note  that  Z(t)  is  a  continuous- time  Markov  process  with  two 
possible  states,  0  and  1.  This  process  is  completely 
defined  under  the  following  assumptions  (Chatfield,  1984); 

Pr[Z(t+At)  -  l|z(t)  »  0]  -  crAt  +  o(At)  (3) 

and 

Pr(Z(t+At)  -  o|z(t)  -  1]  -  /BAt  +  o(At)  (4) 

where  or  is  the  intensity  of  transition  of  a  bubble  from  the 
outside  to  the  inside  of  the  field,  and  ft  from  the  Inside 
to   the  outside.    The  expression  in  the  left-hand  side  of 
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Eqs.  3  and  4,  namely,  Pr[E   E  ],   denotes   the  conditional 

probability  of  event  E  given  event  I  . 

Let 

Pl(t)  -  Pr[Z(t)  -  1]  (5) 

Then,  it  can  be  shown  that  (see  APPENDIX) 

From  the  characteristics  of  the  probe  employed  in  the 
present  work,  we  see  that  the  wave-like  component,  Y(t), 
and  the  stochastic  component,  Z(t),  of  the  recorded 
signals,  X(t),  are  independent  of  each  other;  however,  they 
occur  simultaneously.   Thus,  we  can  write 

X(t)  -  Y(t)Z(t)  (7) 

Auto-Correlation  Function 

The  auto-correlation  function  of  random  data  describes 
the  general  dependence  of  the  datum  at  one  moment  on  that 
at  another  moment  (see,  e.g.,  Bendat,  and  Piersol,  1971). 
Consider  a  sampled  time-history  record  of  X(t);  an 
estimation  for  the  auto-correlation  between  the  values  of 
X(t)  at  times  t  and  (t  +  t)  is  obtained  by  taking  the 
product  of  the  two  values  and  averaging  over  the 
observation  time,  T.  The  resultant  average  product  will 
approach  the  exact  auto-correlation  function  of  X(t)  as  T 
approaches  infinity,  i.e., 


T 
R   (t)  -   lim  £  J   X(t)X(t+T)dt  (8) 

T— ►•    0 

Based  on  the  assumption  that  the  "ergodic  hypothesis"  is 
valid  for  the  system  under  consideration,  the  auto- 
correlation function  is  equal  to  the  corresponding  ensemble 
averaged  value,  i.e., 

Rxx(T)  "  E[X(t)X(t+T)] 


00       » 

J      J   X1X2Q{X1,X2;T)ax.iax2  (9) 


-co     -00 

where  pfx^x^T)  is  the  joint  probability  density  function 

of  X(t)  and  X(t+r). 

Applying  the  definitions   of   the   auto-correlation 
function,   as  given  by  Eq.  9,  and  periodic  component  of  the 

recorded  signals,  Y(t),  given  by  Eq.  1,  we  have 

2* 
Ryy(T)  -   J   A  sin(2jrfQt+e)  A  sin[2rfQ(t+T)+  e]p(e)de 

«=  J-  cos(2irfQT)  (10) 

Since  Z(t)  is  assumed  to  be  a  stationary  random  process, 
the  auto-correlation  of  Z(t)  can  be  derived  as  (see 
APPENDIX) 

Rzz(t)  -  ye"XT  (ii, 

Since  Y(t)  and  Z(t)  are  independent  of  each  other,  Eq.  7 
indicates  that  the  auto-correlation  function  of  X(t)  can  be 
written  as  (see,  e.g.,  Bendat  and  Piersol,  1971) 
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Rxx(T)  ■  Ryy<T)Rzz(T)  <12> 

Substitution  of  Eqs .  10  and  11  into  Eq.  12  gives 

Rxx(T)  "  *T  e"XT  cos<2,rlfoT>  (13> 

Normalizing  the  auto-correlation  function,  we  have 

Rxx(T)     -XT 

Rxx<T)  -  if^loT  *  e    cofi(2'rV)  (14> 


Power  Spectrum 

The  power  spectral  density  function  of  X(t),   6   (u), 

is  defined  as  (see,  e.g.,  Bendat  and  Piersol,  1971) 

2    f* 
Qxx(o)  "  n       J   Rxx(T>  cos<OT>dT  (15) 

Substituting  Eq.   13  into  this  expression  and  integrating, 
we  have 

°XX<«>  -  ^  i    4    °2+2X2r22    2   2  21  <16> 

300        r     u4  +  2(X2-c2)o2+(X2+c2)2 

where 

c  *  2ntQ  (17) 

At  the  stationary  point  in  G^tw),  the  condition 

q    (o)  .  VXA2(-2o)rQ4+2(X2+c2)o2+(X4-2X2c2-3c41 
"  '      [«4+2(X2-c2,«2+(X2+c2)2]2 

holds.   The  value  of  u   satisfying  this  condition  is 
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,.2     2,1/4  ._   ,.2^  2.1/2.1/2  .,_. 

u  -  (X   +  c  )     [2c-(X  +c  )    ]     "  °i         *19) 

We   find  also   that   the  value  of  the  second  derivative  at 
this  point  is  negative,  i.e., 

G   (u) I      <  0 

Thus,  a  single  maximum  of  G   (u)   occurs  at  o  ,   and   the 

corresponding  value  of  this  maximum  is 

2 

q   (0  )  -  ihh i (20\ 

XX*  1*      Aire     2   2  1/2  l*uJ 
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EXPERIMENTAL 

The  facilities  and  procedure  employed  in  carrying  out 
experiments  and  measurements  are  described  in  this  section. 

Facilities 

A  schematic  diagram  of  the  experimental  facilities  is 
shown  in  Figure  1.  The  fluldized  bed  employed  consisted  of 
a  bed  proper,  a  distributor,  and  a  plenum  column.  The  bed 
proper  and  the  plenum  column  were  fabricated  from 
"Plexiglass"  to  permit  visual  observation.  Their  diameter 
was  0.203m  (8  in.),  and  their  heights  were  0.61m  (24  in.) 
and  0.17m  (6.7  in. ), respectively .  The  distributor  was 
perforated  aluminum  plates,  0.00158m  (1/16  in.)  thick  and 
had  164  holes.  The  layout  of  holes  in  the  distributor  is 
shown  in  Figure  2.  Sand  served  as  the  fluldized  particles, 
and  air  at  30-35 *c  the  fluldlzlng  medium.  The  physical 
properties  of  the  sand  and  the  experimental  conditions  are 
summarized  in  Table  1. 

A  pressure  tap  was  installed  on  the  wall  of  the  bed 
column  immediately  above  the  distributor.  The  inside 
opening  of  the  tap  was  covered  with  a  screen  to  prevent  the 
sand  from  entering  the  tap.  The  outside  opening  of  the  tap 
was  connected  to  one  of  the  two  input  channels  of  a 
differential  pressure  transducer  (Enterprise  Model  CJ3D) , 
which  produced  an  output  voltage  proportional  to  the 
pressure  difference   between   two  channels.   The  remaining 
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channel  was  exposed  to  the  atmosphere,  and  thus,  the 
fluctuations  of  pressure  drop  across  the  entire  bed  were 
measured.  The  working  capacity  of  the  transducer  was  ± 
6.90  lcPa  (*  1  psi).  Signal  processing  was  accomplished 
with  the  aid  of  a  Bascom-Turner  digital  recorder,  an  IBM 
PC,  and  a  main  frame  (NAS  6630)  computer. 

Procedure 

For  each  run  of  the  experiment,  the  pressure 
fluctuations  of  the  bed  were  detected  by  connecting  the  tap 
to  the  pressure  transducer.  The  voltage-time  signal 
(corresponding  to  the  pressure-time  signal)  from  the 
transducer  was  fed  to  the  recorder  at  a  selected  sampling 
rate.  The  sampling  rate  was  20  ms.  A  typical  sample 
consisted  of  2000  points.  This  combination  of  sampling 
rate  and  sample  length  ensured  the  capturing  of  the  full 
spectrum  of  hydrodynamic  signal  for  the  gas-solid  fluldized 
bed  (typically  25  Hz).  A  RS-232  Interface  was  used  to  send 
the  signals  from  recorder  to  the  PC.  The  signals  were  then 
sent  to  the  main  frame  computer  via  a  1200  baud  modem  for 
off-line  time-series  analysis. 
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CALCULATIONS 

All  analysis  of  data  were  carried  out  off-line  in  the 
main  frame  computer  utilizing  the  1982  version  of  the  SAS 
(Statistical  Analysis  System)  Package.  The  statistical 
properties  of  the  pressure  fluctuation  signals  calculated 
were  the  mean,  variance,  auto-correlation  function  and 
power  spectral  density  function. 

He  resorted  to  non-linear  regression  analysis  to  fit 
the  model  in  terms  of  the  auto-correlation  function,  given 
by  Eq.  14,  to  the  experimental  auto-correlation  function  by 
manipulating  the  frequency,  tQ,    and  the  intensity,  X.    The 

normalized  power  spectral  density  function,  6   (w)/6   (w., ) , 

was  obtained  by  dividing  Eq.  16  by  Eq.  20.  This  was  done 
to  eliminate  the  parameter  y.  This  spectral  density 
function  ratio  based  on  the  model  was  plotted  over  the 
spectral  density  function,  obtained  from  the  experimental 
data,  to  examine  goodness  of  the  fit.  The  calculations 
were  repeated  for  six  different  superficial  velocities  of 
the  fluidizing  gas. 
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RESULTS  AND  DISCUSSION 

Typical  signals  recorded  from  the  pressure  transducer 
are  traced  In  Figures  3  and  4,  with  superficial  velocity  of 
the  fluldizlng  gas  as  the  parameter.  It  appears  that  each 
tracing  of  the  recorded  signals  at  a  given  superficial 
velocity  contains  wave-like  signals  with  various 
frequencies.  In  other  words,  each  tracing  of  the  recorded 
signals  is  composed  of  a  random  or  stochastic  component  and 
a  wave-like  or  sine  component.  The  variance  or  the 
magnitude  of  each  signal  increases  with  the  increase  in  the 
superficial  velocity.  The  trend  for  the  amplitude  to 
increase  with  the  air  flow  rate  was  expected  since  the 
bubble  size  became  larger  when  air  flow  rate  was  increased. 
The  corresponding  auto-correlation  functions  of  the 
recorded  pressure-fluctuation  signals  at  different 
superficial  velocities  are  illustrated  in  Figure  5  through 
10.  The  shape  of  the  auto-correlation  function  is 
substantially  influenced  by  the  hydrodynamics  of  the  bed 
which  is  directly  related  to  the  pattern  of  bubbling. 
Thus,  analysis  of  the  auto-correlation  function  of  each 
signal  gives  rise  to  statistical  information  concerning  the 
bubble  behavior. 

Comparison  of  the  Model  to  Experimental  Data 

The  auto-correlation  function  based  on  the  model,   Eq. 
14,   has  been   fitted   to   that  obtained  experimentally  by 
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varying  the  frequency,  f  ,  and  the  intensity,  X,   by  means 

of  a  non-linear  regression  analysis.  These  are  Illustrated 
in  Figures  5  through  10.  The  dashed  lines  in  these  figures 
represent  the  model  equation.  We  see  that  the  greater  the 
superficial  velocity,  the  faster  the  decay  of  the  auto- 
correlation function.  This  implies  that  the  intensity  of 
bubble  motion  increases  with  the  increase  in  the  air  flow 
rate ,  i.e.,  the  random  component  of  the  recorded  signals 
becomes  pronounced  compared  to  the  wave-like  signals. 
Figure  11  Illustrates  the  relationship  between  the 
intensity  of  bubbling,  X,  and  the  superficial  velocity  u; 
it  confirms  that  X  Increases  with  u. 

The  frequency  corresponding  to  the  first  major  peak  of 
the  auto-correlation  function  beyond  the  null  delay  (t-0) 
is  identified  as  the  dominant  frequency,  f  .   The  values  of 

fQ     at  various  superficial   velocities  obtained  from  the 

auto-correlation  functions  are  given  in  Table  2;  note  that 
fo  decreases  as  u  increases.   This  may  be  attributed  to  the 

fact  that  the  size  of  the  bubbles  formed  is  bigger  at  a 
higher  air  flow  rate. 

Power  Spectrum 

The  experimentally  obtained  pressure  fluctuation 
signals  have  been  analyzed  and  the  corresponding  power 
spectral   density   functions   are  plotted  as  solid  lines  in 
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Figures  12  through  17  at  different  superficial  velocities 
of  the  fluidizing  medium.  The  most  obvious  and  yet  most 
important  feature  of  the  power  spectral  density  obtained  in 
all  runs  is  the  presence  of  a  sharp  peak.  The  location  and 
other  characteristics  of  the  peaks  depend  on  the  bed 
parameters.  The  concentration  of  the  power  of  the  signal 
into  a  very  narrow  frequency  band  of  less  than  1  Hz 
confirms  the  presence  of  a  periodic  component  in  the 
pressure  signals  (see,  e.g.,  Bendat  and  Piersol,  1971). 
Essentially,  the  same  frequency  is  always  obtained  from  the 
power  spectrum  and  the  auto-correlation  functions  as  shown 
in  Table  2  for  various  superficial  gas  velocities; 
obviously,  the  auto-correlation  and  power  spectral  density 
functions  are  all  consistent  with  one  another.  This 
clearly  demonstrates  the  presence  of  a  periodic  component 
in  the  pressure  fluctuation  signals. 

The  power  spectrum  provides  more  detailed  information 
than  the  auto-correlation  function  on  the  bubble  movement. 
We  observe  in  Figures  12  through  17  that  as  u  Increases, 
additional  small  peaks  appear  in  the  spectral  density 
plots,  thereby  indicating  a  definite  trend  in  the  power 
spectrum,  which  shifts  from  a  predominately  wave-like 
pattern  to  a  predominant  random  pattern.  This  is 
consistent  with  the  observations  made  from  the  auto 
correlation   function  plots   that   as   u   increases,    X 
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(intensity  of  bubbling)  also  increases,  i.e.  the  random 
component  becomes  predominant. 

The  power  spectral   density   function  based  on  the 
model,   Qxx(°) »   is  obtained   from  the   auto-correlation 

function  R   (t)  according  to  Eq.  15.   In  Figures  12  through 

17,  the  experimentally  determined  functions  are  compared 
with  those  calculated  from  the  model  in  Figure  12  through 
17.   Notice  that  the  value  of  the  peak  frequency  (-  u  /2w) 

of  the  power  spectral  density  function  calculated  from  the 
model  is  in  good  agreement  with  that  of  the  frequency,   f  , 

determined  from  the  experimental  data  (see  also  Table  2). 


Relationship  between  the  Frequency  and  Superficial  Velocity 

As  delineated  in  the  introductory  section,  there  seems 

to  be  both  theoretical  and  experimental  variation  in  the 

results  of   the  various  investigators  regarding  the  effect 

of  the  superficial  velocity  of  the  fluidizing  medium  on  the 

frequency   of   pressure   fluctuations.    Hlby   (1967)   has 

proposed  the  following  approximate   relationships  between 

the  bed  oscillation   frequency,  v,      and  u  for   the  two 

limiting  cases  of  laminar  flow  and  turbulent  flow: 

-0.31 
v   cc  u      for  the  laminar  flow  (21) 

and 

u   «  u   for  the  turbulent  flow  (22) 
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According  to  Hiby,  the  higher  the  amplitude  of  oscillation, 
as  a  result  of  higher  porosity,  the  smaller  the  number  of 
layers  which  can  sustain  joint  oscillations.  Nevertheless, 
his   experimental   results  are  not  in  confirmity  with  these 

relations.   He  has  attributed  this  to  the  fact  that  h/d   in 

P 

his  experiment  was  equal  to  only  10;  a  smaller  effective 
height  could  have  meant  an  Increase  in  the  average  natural 
frequency  of  the  coherent  part  of  the  bed,  which  in  the 
process  would  counteract  the  decrease  expected  from  Bqs.  21 
and  22. 

In  the  present  work  the  probe  was  located  about   2   cm 

from  the  distributor;  this  makes  h/d  approximately  200,  if 

h  is  considered  to  be  the  probe  distance  from  the 
distributor.  In  Hiby's  experiments,  the  capacitive  probe 
was  located  at  the  surface  of  the  bed,  thus  making  the 
probe  distance  from  the  distributor  equal  to  the  static  bed 
height. 

Based  on  Hiby's  theory,  an  effort  has  been  made  in  the 
present  work  to  find  an  empirical  relationship  between  u 
and   the  experimentally  determined  f  .   Figure  18  plots  the 

info  vs  inu*  Note  that  the  plot  obeys  the  following 
relationship; 


.-0.623 
o 


f-  «  *  (23) 
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This  approximate  proportionality  appears  to  fall  within  the 
range  of  the  relationships  proposed  by  Hlby  for  the  laminar 
and  turbulent  flows  (Eqs.  21  and  22). 

Fan  et  al.  (1981)  have  proposed  that  the  formation-  of 
large  bubbles  through  bubble  coalescence  in  the  main  body 
or  middle  portion  of  the  bed,  the  formation  of  small 
bubbles  near  the  distributor  and  the  jet  flow  immediately 
above  the  distributor  affect  the  pressure  fluctuations  in 
the  lower  portion  of  a  deep  fluldized  bed.  The  major 
fluctuations  are  caused  by  the  formation  of  large  bubbles 
and  are  transmitted  downward  from  the  middle  portion  of  the 
bed.  The  small  fluctuations,  caused  by  the  jet  flow  and 
the  formation  of  small  bubbles,  are  transmitted  upward  and 
are  superimposed  on  the  major  fluctuations.   The  bed  used 

in  the  present  work  was  relatively  shallow  (H  /D  -  0.542); 

s 

thus,  bubble  coalescence  might  have  occurred   immediately 

above   the  distributor,  forming  larger  bubbles  close  to  the 

probe.   This  could  have  impelled  f   to  decrease  with  the 

o 

increase   in  u.    The   intensity  of  bubbling,  X,  increased 

with  u,  since  the  enhanced  rate  of  air  flow  increased   the 

rate   of   bubble   formation  near   the  distributor,   thus 

increasing  the  randomness  in  the  bed. 
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CONCLUSIONS 

Pressure  fluctuation  signals  from  a  gas-solids 
fluid! zed  bed  have  been  recorded  and  analyzed 
statistically.  The  auto-correlation  and  power  spectral 
density  functions  of  the  pressure  fluctuations  are  all 
consistent  with  one  another;  this  clearly  demonstrate  the 
presence  of  a  periodic  component  in  the  fluctuations.  The 
amplitude  of  the  recorded  signals  increased  with  an 
increase  in  the  superficial  velocity  of  the  fluidizing 
medium,  u  and  it  has  been  found  to  affect  the  pressure 
fluctuations  significantly.    The  frequency,  f  ,  decreased 

and  the  intensity  of  bubble  motion,  X,  increased  when  the 
superficial  velocity  was  increased.  This  phenomena  has 
been  qualitatively  interpreted. 

The  results  of  the  experiment  indicate  that  the 
recorded  signals  are  composed  of  a  wave-like  component  and 
a  stochastic  component.  A  stochastic  model  has  been 
proposed  by  taking  these  two  components  into  account.  The 
auto-correlation  and  power  spectral  density  functions  of 
the  resultant  model  have  been  fitted  to  those  obtained 
experimentally  by  adjusting  the  model  parameters,  the 
dominant  frequency,  fQ,  and  the  intensity  of  transition  of 

bubbles,  X.  The  auto-correlation  and  power  spectral 
density  functions  of  the  model  are  in  good  agreement  with 
those  obtained  experimentally. 
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APPENDIX.   Derivation  of  Eqs.  6  and  11 

Let  us  consider  a  time -continuous  Markov  process  Z(t) 

having  two  possible  states,  0  and   1.  According   to   the 

assumptions   expressed   in  Eqs.  3  and  4  in  the  text,  we  can 

write  the  probability  balance  for  p.(t)  (defined  in  Eq.   5) 

as 

p^t+At)  -  p1(t)[l-/8At]  +  [l-px(t)JaAt  +  o(At)   (A-l) 


or 


Pjit+At)  -  Pi(t) 


- -  -  (a+/8)Pl(t)  +  a  +  Si||i     (A-2) 

Taking  the  limit  as  At  -►  0,  we  obtain 

^|  pA(t)  +  (a+Dp^t)  -  a  <A-3) 

Solving  this  expression  yields 

which  is  Eq.  6  in  the  text.   If 

"l«°>  -  sS?  • 

the  process   is  said   to  be  second  order  stationary  (see, 

e.g.  Chatfield,  1984). 

The  auto-correlation  function  of  Z(t)  is  defined  as 
Rzz(t)  -  E[Z(t)Z(t+r)]  (A-5) 

We  can  see  Immediately  that  the  product  of  Z(t)  and  Z(t+T) 
must  take  on  only  the  values,  0  and  1.  It  will  be  equal  to 
1   if   a  bubble  is  in  the  field  at  both  times  t  and  (t+r); 


21 


otherwise,  it  will  be  zero.  Therefore,  the  expectation  of 
the  product  to  be  equal  to  the  probability  of  this  event  is 
E[Z(t)Z(t+T)]  =  Pr[both  Z(t)  and  Z(t+T)  are  1] 

■  Pr[Z(t)  -  l]Pr[Z(t+T)  »  l|Z(t)  -  1]   (A-6) 

If  we  assume  Z(t)  to  be  a  time  stationary  process,  we  have 
Pr[Z(t+r)  -  ljZ(t)  «*  1]  -  Pr[Z(T)  «  1JZ(0)  «  1)       (A-7) 

Thus,  Eq.  A-6  becomes 

E[Z(t)Z(t+T)]  «=  Pr[Z(t)  -  l]Pr[Z(T)  «  l|z(0)  -  1] 

"  Pi<t>Pn<T>  <A~8> 

Now  considering  the  process  to  be  second  order  stationary 
Eq.  A-4  reduces  to 

Pl(t)  "  £fi  <A~9> 

Substituting  this  expression  into  Eq.  A-8  yields 

E[Z(t)Z(t+T)]  -  ^  Pll(T)  (A-10) 

The  next  step  will  be  to  calculate  the  transition 
probability  P^C*).    We  utilize  the  Kolmogorov  backward 

differential  equations  for  a  two-state  Markov  process  with 
transition  probabilities  or  and  /S,  defined  by  Eqs.  3  and  4 
in  the  text,  respectively.  This  gives  rise  to  the 
following  set  of  differential  equations; 

ft  P01(t)  "  «lPn<*>  "  P0l(t)1  (A"11* 

and 
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|^  Pijtt)  -  /8[P()1(t)  -  P^it)]  (A-12) 

Multiplying  Eq.  A- 11  by  fi   and  Eq.  A-12  by  or  and  adding  the 
resultant  equations  yield 

/8pQ1(t)  +  crp^ft)  -  0  (A-13) 

By  integrating  and  using  the  Initial  conditions 

P01(0)  -  0  and  PjjtO)  -  1,  (A-14) 

we  obtain 

/8p01(t)  -  cr[l-Pll(t)]  (A-15) 

Substituting  Eq.  A-15  into  Eq.  A-12  leads  to 

Pn(t)  -  <or+/B)  [£fi-   P^tt)]  (A-16) 

By  solving   this  expression,   subject   to   the   initial 
condition  p  ,,(0)  -  1,  we  have 

and,  therefore, 

PUm    a+/8  +  a+fi   e  (A-18) 

Substituting  this  expression  in  Eq.  A-10  gives 

.,„tm«rl,  -  ^_  [  _£_  +  _|_  e-<«+*.r  j      (A.19) 

or 

RZZ<T>  -  ^-2  +  -^-2    e'(a+/,>T  (A-20, 

For  the  steady-state  operation,  the  probability  that  a 
bubble  will  enter  the  field  is  identical  to   that   for   the 
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reverse.  Since  the  volume  of  the  field  compared  with  the 
total  volume  of  the  fluidized  bed  is  very  small,  the 
intensity  of  exit  from  the  field,  ft,  is  much  larger  than 
the  intensity  of  transition  from  the  outside  to  the  inside 
of  the  field,  a.  For  a  sufficiently  small  value  of  t,  the 
first  term  on  the  right-hand  side  in  Eq.  A-20  is 
negligible.   Consequently,  Eq.  A-20  may  be  rewritten  as 

R„(T)  -  -^-  e"<a+*)T  -  re"XT  (A-21) 

where 

X  =  <cr+/B) 
and 

aft 


(<*+fi)2 


This  is  Eq.  11  in  the  text. 
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NOMENCLATURE 

A  -  amplitude  of  the  sine  wave  component 

c  *  constant  ■  2irf 


d  =  average  particle  diameter,  ■ 

D  ■  diameter  of  the  bed,  m 

f  ■   cyclic  frequency  of  the  sine  wave  component, 
Hz  or  1/s 

Q   (o)  ■  power  spectral  density  function  of  X(t) 

mtmm 

G   (o. )  ■  maximum  value  of  G„_(«) 

H  «  static  bed  height,  m 

S 

H  =  bed  height  at  minimum  fluidization,  m 
mv 

h  =  height  of  the  probe  from  the  distributor,  m 

PrtEjjEg]  «  conditional  probability  of   event  E±     given 
event  E„ 

p(0)  ■  probability  density  function  of  © 

Pjft)  -  Pr[Z(t)-l],  absolute  probability 

P1;J(t)  -  Pr[Z(t)-j|Z(0)-l],  transition  probability 

P(x1,x2;t)  «   joint   probability  density  function  of  X(t) 
and  X(t+T) 

Rxx*T*  "  aut°-correlation  function  of  X(t) 

* 

Rxx*T*  *  normalized  auto-correlation  function  of  X(€) 

R   (t)  -  auto-correlation  function  of  Y(t) 


Rzz(t)      -  auto-correlation  function  of  Z(t) 
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T  ■  duration  of  the  sample  signal,  s 

t  *  time,  s 

u  =  superficial   velocity   of   the  fluidizing 

medium,  m/s 

u  -  minimum  fluidizing  velocity,  m/s 

nr 

X(t)  -  signal  recorded  from  the  pressure  probe 
Y(t)  ■  sine  wave  component  or  wave-like  signal 
Z(t)        -  stochastic  or  random  component 

Greek  Letters 

a  ■  intensity  of  transition  of  a  bubble  from  the 

outside  of  the  field  to  the  inside 

fi  ■  intensity  of  transition  of  a  bubble  from  the 

inside  of  the  field  to  the  outside 

€  _         =  voidage  at  minimum  fluidization 

mr 

2 

y  «  constant  -  cr£/<cr+/8) 

v  =  oscillation   frequency  of   the   bed    (from 

Hiby's,  1967,  paper),  Hz  or  l/s 

6  ■=  initial  phase  angle  with  respect  to  the   time 

origin,  radians 

r  -  time  shift  variable,  s 

X  =  constant  ■  (a+fi) ,  l/s 

u  ■  angular  frequency,  radian/s 

o  ■  stationary  point  of  G    (u) ,  radian/sec 
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Table  1.   Physical  Properties  of  Sand  and  Test  Ranges  of 
Experimental  Variables. 


Experimental  Variables 


Test  Range 


Sand  Size 
Density  of  Sand 

dp 

u 


■f 


■f 


af 


Sieve  No.  -  3O-d"40 

2620  kg/a3 

0.000491a 
0.25  a/s 

0.54  -  0.98  a/s 
0.47 

0.08-0.09  a 

0.11-0.12  a 
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I    Rotameter 

2  Pressure    gauge 

3  Temperature   gauge 

4  Plenum 

5  Distributor 

6  Bed 
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Figure    1.    Experimental    setup. 
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Figure  3.  Pressure  fluctuation  signal  at  various 
superficial  velocity  of  the  fluidizing 
medium. 
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Figure  4.  Pressure  fluctuation  signal  at  various 
superficial  velocity  of  the  fluidizing 
medium. 
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Figure  5.  Comparison  of  the  auto-correlation  functions 
based  on  the  model  with  that   obtained  from 
the  experimental  data  at  u  r  0.5^32  m/s. 
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Figure  6.  Comparison  of  the  auto-correlation  function 
based  on  the  model  with  that  obtained  from 
the  experimental  data  at  u  =  0.6084  m/s. 
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Figure  7.  Comparison  of  the  auto-correlation  function 
based  on  the  model  with  that  obtained  from 
the  experimental  data  at  u  =  0.65182  m/s. 


37 


DC 


cr 


2.0 


1.5  - 


0.5- 


0.0 


X 

X       - 


0.5 


1.0- 


1.5- 


"2.0 


Experimental    data 

Calculated    curve 

(using    Eq.    I  4) 


'^-^W^ 


1 1 1 1 1 1 1 1 r 

0.0      0.1         0.2       0.3        0.4      0.5        0.6        0.7       0.8        0.9 


r(s) 


Figure  8.  Comparison  of  the  auto-correlation  function 
based  on  the  model  with  that  obtained  from 
the  experimental  data  at  u  =  0.760M  m/s. 
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Figure  9.  Comparison  of  the  auto-correlation  function 
based  on  the  model  with  that  obtained  from 
the  experimental  data  at  u  =  0.8691  m/s. 
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Figure  10.  Comparison  of  the  auto-correlation  function 
based  on  the  model  with  that  obtained  from 
the  experimental  data  at  u  =  0.9777  m/s. 
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Figure  11.  Plot  of  A.  ,  intensity  of  bubbling  vs  u, 
superficial  velocity. 
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Figure  12.  Comparison  of  the  power  spectral  density 
funtion  based  on  the  model  with  that 
obtained  from  the  experimental  data 
at  u  =  0.5432  m/s. 
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Figure  13.  Comparison  of  the  power  spectral  density 
funtion  based  on  the  model  with  that 
obtained  from  the  experimental  data 
at  u  =  0.6084  m/s. 
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Figure  14.  Comparison  of  the  power  spectral  density 
funtion  based  on  the  model  with  that 
obtained  from  the  experimental  data 
at  u  =  0.6518  m/s. 
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Figure  15.  Comparison  of  the  power  spectral  density 
funtion  based  on  the  model  with  that 
obtained  from  the  experimental  data 
at  u  :  0.7604  m/s. 
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Figure  16.  Comparison  of  the  power  spectral  density 
funtion  based  on  the  model  with  that 
obtained  from  the  experimental  data 
at  u  =  0.8691  m/s. 
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Figure    17.    Comparison   of    the   power   spectral    density 
funtion   based   on   the   model   with   that 
obtained    from   the   experimental    data 
at   u   =    0.9777   m/s. 
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Figure  18.  Plot  of  log(f  )  vs  log(u) 
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ABSTRACT 

The  characteristics  of  a  gas-solid  fluidized  bed 
are  influenced  by  complicated  and  stochastic  phenomena, 
e.g.,  jetting  and  bubbling  of  the  fluidizing  medium  and 
the  motion  of  the  fluidized  particles.  In  this  study 
pressure  fluctuations  due  to  the  bubble  movement  were 
measured  by  means  of  a  pressure  probe  over  a  range  of 
time.  The  resultant  time-series  have  been  analyzed  by 
determining  their  auto-correlation  function  and  power 
spectrum  with  emphasis  on  the  examination  of  the  effect 
of  the  superficial  velocity  of  the  fluidizing  gas.  A 
stochastic  model  of  bubble  motion  in  a  fluidized  bed 
has  been  developed.  This  model  visualizes  the  bubble 
motion  in  a  fluidized  bed  to  consist  of  the  random 
movement,  generating  irregular  signals,  and  the  linear 
movement,  generating  wave-like  signals.  A  theoretical 
auto-correlation  function  and  a  power  spectral  density 
function  have  been  derived  based  on  the  model. 


